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filter  processor,  the  fast  Fourier  transform  (FFT),  and  MTIs  cascaded  with  the  FFT.  The 
processors  are  compared  on  the  basis  of  their  clutter  rejection,  signal  processing  losses,  and 
minimum  detectable  velocity  (MDV)  capabilities. 
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A  MATCHED  FILTER  DOPPLER  PROCESSOR 
FOR  AIRBORNE  RADAR 


INTRODUCTION 

Two  complications  inherent  in  the  detection  of  aircraft  from  an  airborne  moving  platform  are  the 
level  of  clutter  rejection  required  (70  dB  or  better  is  the  goal  in  this  research)  and  the  difficulties  arising 
from  the  platform  motion  which  both  shifts  the  clutter  spectrum  and  increases  its  width. 

Another  problem  with  Doppler  frequency  processing  is  the  minimum  detectable  velocity  (MDV), 
which  introduces  blind  regions  in  the  surveillance  region.  The  width  of  the  blind  region  is  directly 
related  to  the  MDV  (i.e.,  the  higher  the  MDV  the  wider  the  blind  region).  Because  of  such  blind 
regions,  there  exists  the  possibility  of  a  target  penetrating  a  surveillance  region  by  flying  a  course  that  is 
nearly  tangential  to  the  radar  platform. 

In  the  research  reported  here,  it  is  shown  that  a  matched  Alter  signal  process  or  [1]  will  produce 
the  necessary  level  of  clutter  rejection  to  distinguish  moving  from  Axed  targets.  The  results  are  com¬ 
pared  to  those  obtained  using  a  fast  Fourier  transform  (FFT)  processor  and  with  those  obtained  using  a 
moving  target  indicator  (MTI)  cascaded  with  an  FFT  Doppler  Alter  bank  [2]. 

THE  MATCHED  FILTER  CONCEPT 

Applebaum's  adaptive-antenna  algorithm  13]  can  be  applied  to  the  control  of  a  bank  of  Doppler 
Alters  [4],  From  [3]  and  [4]  the  optimum  Alter  weight  vector  for  a  matched  Alter  (i.e.,  when  the  signal 
is  known  completely  so  that  the  signal  vector  5  can  be  deAned)  is 

-  Mr'S*.  (1) 

where  A//-1  represents  the  inverse  interference  covariance  matrix  and  the  superscript  (*)  denotes  the 
complex  conjugate. 

In  [1]  a  generalized  Alter  is  derived  by  assuming  that  the  signal’s  Doppler  shift  is  not  precisely 
known  but  is  known  only  within  some  interval.  The  entire  Doppler  region  of  interest  is  covered  with  a 
bank  of  Alters,  one  for  each  interval. 

For  very  small  pulse  widths  compared  to  the  time  between  pulses  [S],  the  received  signal  has  the 

form 

s(r)  -  8 (/  -  kT)  exp  (Jlnfjt),  (2) 

where  Tis  the  time  between  transmitted  pulses, 

{0,  t  *  kT 
1  t  -  kT. 

Assuming  the  Doppler  shift  ft  has  a  uniform  probability  density  function  over  an  interval  of 
length  B  centered  at  fh  then 
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_L  ,  B 

B‘  Ji~  j  ^  Jd  <  Ji+  2 

f(/rf)  “  o  elsewhere,  ^ 

where  ft  is  the  center  of  the  ith  filter  and  is  given  by  /,  ”  0  —  D/NT,  i  —  1,  2 .  Af  and  5  — 

1/ AT  where  Ar  is  the  number  of  pulses  integrated. 

Figure  1  (using  an  8-pulse  processor  as  an  example)  shows  the  principal  difference  between  the 
Applebaum  algorithm  and  the  design  presented  in  (1)  and  applied  here  to  an  airborne  system.  In  Fig. 
la  each  filter  is  optimized  (i.e.,  matched)  to  one  particular  Doppler  frequency.  In  Fig.  lb,  the  design 
used  here,  the  filter  is  optimized  to  detect  moving  targets  that  have  Doppler  frequencies  that  fall  any¬ 
where  within  one  of  the  intervals  (i.e.,  it  is  matched  to  the  expected  value).  Thus,  f  represents  the 
center  of  the  /th  filter  with  the  first  filter  centered  on  zero  Doppler  as  shown  in  Fig.  1. 
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Fig.  1  —  Basis  for  the  design  of  a  bank  of  Alters  to  reject  the  clutter 
and  detect  a  Urget  with  (a)  a  Doppler  shift  at  one  of  the  equally 
spaced  points  or  with  (b)  a  Doppler  shift  that  may  be  anywhere 
within  one  of  the  equal  intervals  of  the  MTI  detection  region 


The  expectation  of  the  signal  with  respect  to  the  unknown  Doppler  shift  is 

7(f)  -  £ls(f)}  -  80  -  kT)  £  P<Ji)  exp  02vf4t)df4. 

The  integral  above  is  the  inverse  Fourier  transform  of  the  probability  density  function  of  the  Doppler 
shift  which  yields 

JO)  -  80  -  kT)  exp  (J2wf,t)  "'"J".  (4) 

All  the  parameters  of  Eq.  (4)  are  known  except  k. 

An  N- pulse  Doppler  filter  requires  N  samples  of  7(f).  We  must  choose  N  sequential  samples 
from  the  infinite  set  described  by  Eq.  (4).  The  filter  response  to  7(f)  is  maximized  if  the  samples 
selected  are  centered  around  the  peak  of  the  (sin  x)/x  curve  of  Eq.  (4). 
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The  /th  component  of  the  signal  vector  for  the  i th  filter  is 

2wU  -  1) 


s,  -  exp  \j 


N 


I  - 


N  +  1 


v_  N  +  1 
N  2 


(5) 


for  /-  1,  2 . N  and  /' -  1,  2 . N. 


Since  Eq.  (5)  represents  a  completely  known  signal  vector,  it  can  be  used  in  the  Applebaum  algo¬ 
rithm 

-  Mr' S*  (6) 

where  the  /th  component  of  S*  is  given  by  the  complex  conjugate  of  Eq.  (S).  By  using  the  expected 
value  of  the  signal  rather  than  the  known  value,  we  do  not  have  a  matched  filter  in  the  strict  meaning 
of  the  phrase.  What  we  actually  have  is  a  filter  matched  to  the  expected  value  of  the  unknown  signal 
vector.  Therefore,  weight  vectors  can  be  computed  by  using  Eq.  (6). 

One  then  determines  the  platform  motion  spectrum.  The  platform  motion  spectrum  (Fig.  2)  con¬ 
sists  of  samples  that  are  computed  from  the  gain  of  the  antenna  pattern  at  a  given  angle  relative  to 
boresight  multiplied  by  the  backscatter  of  the  clutter  taken  in  equal  steps  of  frequency  (not  angle)  and 
converted  to  the  azimuth  angle. 


For  a  0°  dive  angle,  applying  Eq.  (A3)  gives, 

,  2v, 

fd  -  -jjt-  cos  <t>  cos  a 

where  <f>  is  the  depression  angle  of  the  antenna  relative  to  the  velocity  (v.)  and  a  is  the  azimuth. 


(7) 


Of  it 


-20 


DOPPLBt  SHIFT/ RADAR  FRF 
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Fig.  2  —  Clutter  spectrum  due  to  platform  motion;  S  nmi  (30  000  ft)  altitude, 
30-kHz  pulse  repetition  frequency,  azimuth  is  90* 
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By  Eq.  (7), 

2v 

/,  +  i»A/rf  -  — ^  cos  ♦  cos  (aa  +  a,)  (8) 

where  afl  is  the  azimuth  of  boresight  relative  to  the  velocity,  f0  is  the  Doppler  frequency  on  boresight, 
n  is  an  integer  representing  the  number  of  equal  frequency  steps  with  respect  to  boresight,  and  a„  is 
the  azimuth  (relative  to  boresight)  corresponding  to  the  frequency  represented  by  the  left-hand  side. 


Solving  Eq.  (8)  for  a,  we  get 

a* 


If  we  let  k 


A 

2v,  cos  4  ’ 


then 


cos  1 


( fa  +  nA/,)A 

2v,  cos  4> 


Ota. 


cos~1((/«  +  hA/j)k!  -  a„, 


(9) 


where  n  *  ^mm,  *  .  2,  1,  0,  1,  2,  ... ,  Rmw 

Choosing  a  A  fd  to  assure  an  adequate  number  of  points,  we  determine  both  nmin  such  that 

2v 

/«  -  »mi COS  (10) 

and  n„ „x  such  that 

/.  +  cos  +.  (11) 

/rmin  &  nmx  except  where  a.  —  90“ . 


The  toUl  number  of  points  for  a  given  4  and  a„  is 

II  ■*  Amin  “I"  nmix  ^  !•  (12) 

Thus,  we  can  determine  the  azimuth  angles  with  respect  to  boresight  by  Eq.  (9),  e.g., 

Ot min  “  COS  [(/j  *  tlnjgA/^)l(]  ® I 


a_j  —  cos-1 1(/,  -  2A/rf)ic]  -  a, 

a- 1  -  COS'1 1 (/,  -  1A/,)k]  -  a, 

a0  “  cos~‘[(/,  +  0A/*)k)  -  a,  (boresight) 
<*i  —  cos'1 ((/,  +  1A/*)*J  -  a, 
a2  —  cos'1  [(/,  +  2A  /rf)#tl  —  a. 


<*««  -  cos"1 1(/«  +  ~  a. 
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where  «min  represents  the  azimuth  angle  (relative  to  boresight)  where  the  minimum  Doppler  frequency 
>  — 2vp/k  cos  <j>  occurs,  0^,  represents  the  azimuth  angle  (relative  to  boresight  *nere  the  maximum 
Doppler  frequency  <  2 vp/x  cos  <l>  occurs,  and  the  other  subscripts  designate  the  azimuth  angles  (rela¬ 
tive  to  boresight)  corresponding  to  the  equivalent  number  of  frequency  steps  with  respect  to  boresight 
(0  for  boresight,  negative  for  frequencies  less  than  boresight,  and  positive  for  frequencies  greater  than 
boresight). 

Next,  to  determine  the  area  of  the  clutter  for  each  sample,  apply 

Ae  —  -y-  sec  2R  tan  1^‘JJ’  03) 

where  R  is  the  slant  range  corresponding  to  the  depression  angle  </>,  is  the  grazing  angle,  c  is  the  pro¬ 
pagation  speed,  t  is  the  pulse  width,  and  40  is  a  function  of  a. 

40  is  computed  by  one-half  the  sum  of  the  two  angular  widths  (angular  widths  are  determined  by 
subtracting  each  consecutive  pair  of  results  from  Eq.  (9))  adjacent  to  an  azimuth  angle  computed  by 
Eq.  (9). 

Finally, 


A  similar  procedure  is  performed  in  the  Appendix  for  filter  widths  and  the  area  of  the  clutter  for 
the  footprint.  The  variation  of  azimuth  angular  widths  when  equal  frequency  steps  are  used  can  be 
seen  clearly  in  Fig.  (A2).  However,  the  frequency  steps  used  in  the  Appendix  are  not  the  same  as 
those  applied  here.  Here,  we  require  a  relatively  large  number  of  points  (i.e.,  much  smaller  4 fds). 

The  interference  spectral  density  due  to  internal  motion  (the  clutter  has  a  Gaussian  frequency 
spectrum;  the  noise  has  a  "white"  frequency  spectrum;  the  clutter  and  noise  both  have  a  Gaussian 
amplitude  probability  density  function;  normalized  to  unity  Ec  +  Ew  -  1)  is  given  by 

*  — _ I  Wt“Mc)2l  .  -  -[..//X  ..I,  ill  / 1 


P,(Jj)  "  V2rr 


exp  I- 


+  EWT  «(/rf)-u/a-y 


where  the  mean  nc  implies  relative  motion  between  the  radar  platform  and  the  clutter,  <rc  (the  standard 
deviation  of  the  clutter  spectrum)  is  a  measure  of  the  bandwidth  of  the  clutter  spectrum,  T  is  the  time 
between  transmitted  radar  pulses,  Ew  is  the  power  of  the  "white"  spectrum  Ec  is  the  power  of  the  Gaus¬ 
sian  spectrum,  and 

{0.  /,  <  0 
m{fJ  “  |l.  /,  >  0. 


The  internal  motion  spectrum  is  computed  by  Eq.  (15)  by  using  the  same  frequency  steps  as  one 
does  in  computing  the  platform  motion  spectrum. 

The  interference  covariance  matrix,  Mf,  may  be  determined  by  first  convolving  the  spectrum  of 
the  internal  motion  plus  noise  (Eq.  15)  with  the  platform  motion  spectrum  (Fig.  2)  to  yield  the  total 
spectrum. 

Second,  take  the  inverse  discrete  Fourier  transform  of  the  total  spectrum  [51.  This  produces  a 
discrete  autocorrelation  of  the  delay  (r)  where  r  is  defined  by 

r  -  nT,  - N  <  n  <  N  and  n  €  (integers).  (16) 
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Finally,  the  interference  covariance  matrix  elements  may  be  determined  as  presented  in  the 
Appendix  of  (5). 

Matched  FUter 

Figure  3  shows  the  shape  of  filter  number  33  of  the  matched  filter  processor  with  optimal  weights 
defined  by  Eq.  (6).  The  internal  motion  is  small  compared  to  the  platform  motion.  Therefore,  the 
filter  is  approximately  matched  to  the  platform  motion  spectrum  of  Fig.  2.  The  reciprocal  relationship 
between  the  filter  response  of  Fig.  3  and  the  spectrum  of  Fig.  2  can  be  seen.  The  loss  is  about  0.1  dB. 
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Fig.  3  —  Filter  No.  33  of  t  64-putee  matched  Biter  processor  with  optimal 
weights.  This  Biter  is  optimized  to  the  'input*  spectrum  of  Fig.  2. 


Figure  4  shows  the  improvement  factor  (IF)  of  the  matched  filter  processor  for  this  example. 
The  clutter-to-noise  ratio  is  6S  dB.  The  MDV  is  the  velocity  at  which  we  can  achieve  a  70-dB  improve¬ 
ment  factor  (explained  later).  For  this  processor  it  is  about  80  knots. 

Cascaded  MTI  aad  Matched  Filter 

For  a  cascaded  MTI-Matched  filter  configuration,  the  filter  is  matched  to  the  covariance  matrix 
corresponding  to  the  output  spectrum  shown  in  Fig.  5  which  shows  the  results  of  a  single  MTI  canceller 
(main  lobe  clutter  is  reduced  by  about  34  dB). 

Figure  6  shows  the  IF  for  this  processor  with  optimal  weights  defined  by  Eq.  (6).  The  clutter-to- 
noise  ratio  is  65  dB.  The  MDV  is  about  80  knots  and  the  loss  is  about  0.1  dB.  The  IF  is  about  5  dB 
teas  than  shown  in  Fig.  4. 
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Fig.  6  —  Improvement  Factor  for  a  64-pulse  Doppler  processor  using  the  matched  filter 
algorithm  preceded  by  a  single-canceller  MTI.  The  clutter-to-noise  ratio  is  6S  dB. 


Since  the  matched  filter  without  the  MTI  is  a  much  simpler  processor  design  and  since  there  is  no 
degradation  of  performance,  it  appears  to  be  the  better  configuration  However,  an  MTI-matched  filter 
system  using  a  single-canceller  requires  about  six  bits  less  precision  in  the  processor  as  shown  by  its 
cancellation  ratio  (CR  on  Fig.  5)  of  about  35  dB.  This  may  be  desirable  for  some  applications  because 
it  would  require  less  total  hardware. 

Cascaded  DPCA  Canceller  and  Matched  Filter 

Two  commonly  used  methods  of  correcting  for  platform  motion  are  time  averaged  clutter 
coherent  airborne  radar  (TACCAR)  and  displaced  phase  center  antenna  (DPCA).  TACCAR  .16] 
removes  the  average  Doppler  shift  due  to  platform  motion,  and  DPCA  [7]  corrects  for  the  spectral 
spreading  also  due  to  the  platform  motion.  Figure  7  shows  the  results  when  a  DPCA  canceller  is  cas¬ 
caded  with  a  64-pulse  matched  filter  processor.  The  clutter-to-noise  ratio  is  65  dB.  The  IF  is  about  the 
same  as  that  shown  in  Fig.  4  for  the  matched  filter  alone.  However,  the  MDV  is  about  45  knots  vs 
about  80  knots.  This  is  a  remarkable  improvement  in  the  MDV.  Although  this  processor  requires  less 
processing  hardware  for  the  same  reason  as  the  MTI,  it  introduces  antenna  complications  in  that  usually 
three  antenna  patterns  must  be  generated  [7], 

COMPARISON  WITH  OTHER  TECHNIQUES 

Current  radar  Doppler  processing  techniques  to  distinguish  moving  from  fixed  targets  consist  of: 

1.  The  moving  target  indicator  (MTI), 

2.  The  weighted  fast  Fourier  transform  (FFT),  and 

3.  The  cascaded  MTI/FFT. 
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Fig.  7  —  Improvement  factor  for  a  DPCA  canceller  cascaded  with  a  64-pulse 
matched  filter  processor.  The  clutter-to-noise  ratio  is  65  dB. 


Criteria  used  to  evaluate  the  performance  of  a  particular  Doppler  processor  are:  the  improvement 
factor  (IF),  which  is  a  measure  of  the  clutter  rejection  capability;  the  signal  processing  loss,  which  is  a 
measure  of  the  loss  in  signal  to  noise  ratio  from  the  ideal;  and  the  minimum  detectable  velocity 
(MDV),  which  is  a  measure  of  the  ability  to  detect  low  velocity  targets  or  targets  whose  velocity  vector 
is  nearly  tangential  to  the  direction  of  the  radar. 

The  limitations  and  disadvantages  of  the  former  processors  are  primarily  in  the  areas  of  the  signal 
processing  loss  and  MDV  response. 

Feedback  MTI 

The  conventional  MTI  canceller  delays  the  returns  of  a  given  transmitted  pulse  and  subtracts 
them  from  the  returns  of  the  next  transmitted  pulse.  N  cascaded  cancellers  are  equivalent  to  an 
(N  +  1)  sample  transversal  filter  with  weights  corresponding  to  the  Mh-degree  bionomial  coefficients 
with  alternating  signs  [5]. 


In  general,  the  binomial  weights  are 


where 


(-I)'  i  '  1 . Nc 


i  J  “  /((jv  _  1) i  anc*  is  the  number  of  cancellers. 


For  a  triple  canceller  ( N  -  3),  the  binomial  weights  would  be  (1,  -3,  3,  -1),  which  when  nor¬ 
malized  by  dividing  by  the  absolute  value  of  the  largest  binomial  coefficient  would  be  (0.333,  —1.0,  1.0, 
-0.333). 
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Figure  8  shows  the  filter  gain  relative  to  the  maximum  fitter  gain  vs  Doppler  frequency  relative  to 
the  radar  pulse  repetition  frequency  (PRF)  for  1  to  7  MTI  cancellers  [2]. 

By  cascading  MTI  cancellers  the  clutter-rejection  null  (centered  around  0  and  1  on  the  horizontal 
axis)  is  broadened  as  shown  by  Fig.  8.  For  a  given  clutter  spectral  distribution,  the  broader  null  yields 
a  better  IF  as  shown  by  Fig.  9.  The  penalty  paid  for  the  gain  in  the  IF  is  a  poor  MDV  response  [8]  as 
welt  as  wide  blind  velocity  regions  that  occur  when  the  Doppler  frequency  is  a  multiple  of  the  PRF. 
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Fig.  8  —  Normalized  power  transfer  function  for  an  MTI  having  the 
indicated  number  of  cancellers 
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Another  concurrent  problem  with  the  poorer  MDV  response  is  the  increase  of  the  loss  (Fig.  10) 
as  the  number  of  cancellers  cascaded  increases  where 

Loss  -  10  log  ( Ne  +  1).  (18) 

The  loss  corresponds  to  the  loss  of  an  MTI.  However,  if  the  MT1  is  followed  by  coherent  or  non¬ 
coherent  integrators,  this  loss  can  be  reduced,  but  not  eliminated. 

One  method  of  improving  the  poor  MDV  response  of  the  MTI  is  by  shaping  the  skirl  of  the  filter 
by  various  weighting  schemes  as  shown,  for  example,  in  Fig.  11  from  [9]  where  1/Tis  equal  to  the 
PRF.  This  method  improves  the  MDV  response  but  the  shaping  requires  longer  processing  times 
(which  implies  more  pulses  processed)  which  in  turn  increases  the  loss. 


Fig.  10  —  toss  for  an  MTI  having  the  indicated  number  of  cancellers 


FREQUENCY 

Fig.  1 1  —  Amplitude  responses  for  three  MTI  delay-line  cancell¬ 
ers,  (1)  Classical  3-pulse  canceller,  (2)  5-pulse  delay-line  can¬ 
celler  with  ’optimum'  weights,  and  (3)  15-pulse  Chebyshev 
design.  (After  Houts  and  Burlage  (9).) 
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Coherent  Integrator  (FFT) 

The  weighted  FFT  gives  a  better  MDV  response  than  the  MTI  [2,  8),  but  it  also  produces  2  to 
3-dB  losses  because  of  the  weighting  function  which  lowers  the  filters’  sidelobes  everywhere,  not  just  in 
the  vicinity  of  the  clutter  (see  Fig.  12). 
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Fig.  12  —  Filler  No.  33  of  a  64-pulse  FFT  with  -70  dB 
Chebyshev-weights  —  loss  ~  2.2  dB 


Heavy  weighting  (to  control  the  level  of  the  sidelobes)  is  needed  for  the  required  IF.  The  result 
is  a  wider  main  lobe  with  its  concurrent  poorer  MDV  response  because  the  filters  closer  to  the  main 
lobe  clutter  reject  less  clutter  causing  the  output  clutter  residue  to  be  dominated  by  the  main  lobe 
clutter.  The  effect  of  this  domination  is  a  reduced  detectability  for  the  low  velocity  targets  in  the  filters 
close  to  the  main  lobe  clutter  [8]. 

For  comparison  with  the  matched  filter  processors  formerly  presented.  Fig.  12  shows  filter  number 
33  of  a  64-pulse  FFT  with  -70  dB  Chebyshev  weights  (chosen  as  a  result  of  the  required  IF  of  70  dB). 
The  loss  is  about  2.2  dB. 

The  IF  for  a  64-pulse  FFT  with  a  clutter-to-noise  ratio  of  65  dB  and  -70-dB  sidelobes  is  shown  in 
Fig.  13.  The  IF  for  filter  number  33  is  about  75  dB;  the  MDV  is  about  120  knots. 

Figure  14  shows  the  input  and  output  spectra  of  filter  number  33  for  the  FFT  64-pulse  processor. 
The  figure  shows  that  the  output  spectrum  is  still  dominated  by  the  main  lobe  clutter  residue.  We  can 
further  reduce  the  main-lobe  clutter  residue  by  preceding  the  FFT  processor  with  one  or  more  moving 
target  indicators  (MTIs). 
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Cascaded  MTI  aad  FFT 

The  performance  of  coherent  integrators  against  white  noise  is  well  established.  Since  the  perfor¬ 
mance  in  colored  noise  (or  clutter)  is  irregular  and  unpredictable  [2],  we  will  partially  prewhiten  the 
clutter  input  to  the  coherent  integrator  by  preceding  it  with  an  MTI  canceller. 

Figure  IS  shows  the  IF  for  a  double-canceller  MTI  cascaded  with  a  64-pulse  FFT  with  uniform 
weights.  With  uniform  weights,  the  loss  is  less  than  0.1  dB.  The  MDV  is  increased  to  about  170  knots 
due  to  the  poor  MTI  low  velocity  response. 


FiS  15  —  Improvement  factor  for  ■  double-canceller  MTI  cascaded  with  a  64-pulse 
coherent  intesralor  with  uniform  weights  (unweighted) 


DESIGN  PROCEDURE 

The  matched  filters  could  be  designed  to  give  the  best  performance  by  computing  the  weights  for 
the  clutter  spectrum  at  each  antenna  pointing  angle.  This  would  probably  need  to  be  implemented 
adaptively.  An  alternative  method  would  be  to  compute  a  finite  number  of  sets  of  weights  at  specific 
angles  (e.g.,  0s,  45s,  and  90”).  The  set  corresponding  to  the  angle  closest  to  the  horizontal  pointing 
angle  could  then  be  picked  and  used.  The  Doppler  spectral  spreading  due  to  range  ambiguities  can  be 
compensated  for  with  quadratic  phase  modulation.  However,  for  this  analysis  we  will  take  only  a  single 
range  cell. 

The  matched  filters  of  the  example  presented  here  were  designed  for  the  widest  clutter  spectrum 
and  highest  dutter-to-noise  ratio  (CNR)  expected.  The  worst  case  weights  (i.e.,  for  a  horizontal  point¬ 
ing  angle  of  90”  which  maximizes  the  effects  of  platform  motion)  were  computed  and  stored.  These 
weights  were  then  used  at  all  other  horizontal  pointing  angles. 
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Figure  16  shows  the  target  outputs  of  an  airborne  radar  model  (Xs)  superimposed  on  the  matched 
filter  improvement  factor  curves  for  the  worst  case  mentioned  above.  The  radial  velocities  of  the  tar¬ 
gets  that  were  modeled  were  selected  to  be  at  the  center  of  filters  2,  3,  4,  5,  7,  9,  17,  and  33.  The  hor¬ 
izontal  pointing  angle  for  both  outputs  was  90s.  The  CNR  for  both  outputs  was  65  dB.  As  shown  in 
Fig.  16,  the  modeled  outputs  match  the  filter  responses  rather  well.  The  circled,  modeled  output 
represents  the  performance  using  the  skirt  of  filter  3  rather  than  the  center  of  filter  2  for  the  target  that 
is  centered  in  filter  2.  Since  the  skirt  of  filter  3  (in  this  case)  outperforms  filters  1  and  2,  a  good  case 
could  be  made  for  discarding  the  results  of  filters  1  and  2  (and  also  filter  64). 
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Fig.  16  —  Matched  Alter  modeled  target  outputs  (Xs)  superimposed  on  the  matched  Alter 
improvement  factor  curves.  The  horizontal  pointing  angle  is  90*  (for  both  targets  and  Alter 
design).  The  clutter-to-noise  ratio  is  65  dB  (for  both  targets  and  Alter  design). 


Figure  17  shows  the  modeled  outputs  with  a  55-dB  CNR.  Both  outputs  still  represent  a  90*  hor¬ 
izontal  pointing  angle.  The  IF  is  about  10  dB  less.  However,  the  IF  for  filter  2  is  about  the  same  (i.e., 
about  60  dB).  With  a  CNR  of  55  dB  we  require  10  dB  less  cancellation  (i.e.,  60  dB  rather  than  70  dB). 
We  can  conclude  that  our  matched  filter  processor  using  the  90*  weights  and  designed  for  a  65-dB  CNR 
performs  at  least  as  well  in  terms  of  the  MDV  against  the  55-dB  CNR  clutter  as  it  did  against  the  65-dB 
CNR  clutter. 
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Fig.  17  —  Modeled  target  output  (Xs)  superimposed  on  the  improvement  factor  curves.  The 
horizontal  pointing  angle  is  90°  (for  both  targets  and  Alter  design).  The  clutter-lo-noisc  ratios 
are  SS  dB  (targets)  and  65  dB  (Alter  design). 


Figure  18  shows  the  modeled  outputs  with  a  CNR  of  55  dB  and  a  horizontal  pointing  angle  of  0” 
superimposed  on  the  filter  response  matched  to  a  CNR  of  6S  dB  and  horizontal  pointing  angle  of  90°. 
The  IF  here  is  essentially  the  same  as  that  in  Fig.  1 7  except  for  the  filters  near  the  clutter,  where  it  is 
actually  better  if  the  skirt  of  filter  3  is  used  rather  than  the  center  of  filter  2. 

As  shown  by  the  modeled  outputs  in  Fig.  18  where  we  have  10  dB  less  CNR  and  the  90s  weights 
used  at  0°  horizontal  pointing  angle  (the  worst  possible  utilization  of  the  worst  case  weights),  the  per¬ 
formance  in  terms  of  the  MDV  is  just  as  good  and  the  cancellation  level  required  is  still  achieved.  One 
can  conclude  then  that  this  particular  set  of  weights  will  suffice  for  this  example.  Also,  we  have  com¬ 
pared  the  matched  filter  using  these  weights  to  the  FFT  processors  earlier  in  this  report.  We  found  that 
the  matched  filter  processor  outperforms  these  configurations  without  adapting  the  weights  for  each 
antenna  pointing  angle.  Therefore,  the  implementation  of  this  processor  is  not  significantly  more  com¬ 
plex  than  the  FFT  processors. 
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Fig.  18  —  Modeled  target  (Xs)  superimposed  on  the  improvement  factor  curves.  The  hor¬ 
izontal  pointing  angles  are  0°  (targets)  and  90°  (Alter  design).  The  clutter-to-noise  ratios  are 
SS  dB  (targets)  and  6S  dB  (Alter  design). 


SUMMARY  AND  CONCLUSIONS 

We  have  considered  five  processor  configurations  in  this  research  to  be  applied  to  an  airborne 

radar: 

1 .  FFT  processor  with  -70  dB  Chebyshev  weights, 

2.  Double-canceller  MTI  cascaded  with  an  FFT  processor, 

3.  Single-canceller  MTI  cascaded  with  a  matched  filter  processor, 

4.  Matched  filter  processor,  and 

5.  DPCA-canceller  cascaded  with  a  matched  filter  processor. 

Using  a  design  example,  our  goal  was  to  achieve  70  dB  or  better  clutter  rejection  with  minimal 
loss  and  a  minimum  detectable  velocity  (MDV)  of  100  knots  or  less. 
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A  comparison  of  these  five  configurations  yields  the  following  approximate  results: 


Type  of  Processor 

MDV 

(knots) 

Loss 

(dB) 

1.  FFT 

120 

2.2 

2.  MTI-FFT 

170 

0.1 

3.  MT1-MATCHED 

80 

0.1 

4.  MATCHED 

80 

0.1 

5.  DPCA-MATCHED 

45 

0.1 

Our  conclusions  are: 

1.  70-dB  or  better  clutter  cancellation  is  achievable  with  all  configurations  for  our  design 
example. 

2.  Configurations  3,  4,  and  S  each  have  a  low  enough  MDV  to  meet  the  criteria  set  by  our 
design  example. 

3.  Configurations  3  and  4  give  the  same  results.  Configuration  3  requires  less  processing 
hardware  because  the  cancellation  ratio  of  the  MTI-canceller  reduces  the  dynamic  range 
required  for  the  matched  filter.  Configuration  4  is  a  more  straightforward  design  which  is 
based  on  the  spectrum  due  to  platform  motion,  whereas  the  design  of  Configuration  3  is 
based  on  the  platform  motion  spectrum  after  the  MTI  cancellation. 

4.  Configuration  S  produced  the  lowest  MDV  with  overall  results  as  good  as  the  other 
configurations  but  Configuration  5  requires  three  antenna  patterns  to  be  generated  which 
complicates  the  antenna  design. 

5.  In  Configuration  4,  the  filters  were  matched  to  a  wide  clutter  spectrum  (at  broadside)  and 
performed  at  least  as  well  for  a  narrow  clutter  spectrum  (straight-ahead)  in  terms  of  the 
MDV.  If  the  filters  were  perfectly  matched  we  could  achieve  a  better  signal-to-noise  ratio 
(a  reduced  loss),  but  the  difference  between  such  a  reduced  loss  and  the  0.1 -dB  loss  of 
our  design  example  would  be  insignificant  for  most  applications,  especially  when  the 
added  complications  are  considered. 

6.  The  primary  motive  for  implementing  Configuration  4,  the  matched  filter  processor,  is  its 
simplicity  and  the  fact  that  only  one  set  of  weights  needs  to  be  computed  and  applied  at 
all  pointing  angles. 
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Appendix 

AIRBORNE  RADAR  CLUTTER  MODEL 

Clatter  Spectrum 

For  an  airborne  radar,  curves  of  constant  Doppler  shift  (isodops)  can  be  determined  by  the 
Doppler  frequency  equation 

2v 

fd  -  cos  0,  (Al) 

where  fd  represents  the  Doppler  frequency,  is  the  magnitude  of  the  platform  velocity  vector,  A  is 
the  transmitted  wavelength,  and  0  is  the  angle  measured  from  the  velocity  vector  to  the  line  of  sight  of 
an  earth  point  (Fig.  Al).  From  the  coordinate  system  defined  in  Fig.  (Al), 

cos  0  “  cos  D  cos  ^  cos  a  4-  sin  D  sin  <£,  (A2) 

and 

2v 

fd  -  -  (cos  D  cos  <f>  cos  a  +  sin  D  sin  4>).  (A3) 


TO  EARTH'S 
CENTER 

Fif.  A I  —  Geometry  of  in  airborne  radar  showing  the  relationship  of 
azimuth  angle  (a),  depression  angle  (♦),  and  dive  angle  (D) 


In  modeling  the  clutter,  isodops  are  determined  by  Eq.  (A3).  The  total  clutter  level  in  a  given 
Doppler  filter  is  determined  by  summing  the  clutter  in  the  ambiguous  range  cells  along  the  isodop  cen¬ 
tered  on  the  filter.  The  isodops  centered  on  the  filters  are  determined  by  using  frequency  steps  of 

A/d  -  where  T  is  the  interpulse  period  and  N  is  the  number  of  pulses  to  be  processed. 

NT 

Figure  A2  shows  the  isodops  for  64  filters  and  0°  dive  angle.  Outermost  circles  represent  the  hor¬ 
izon  for  a  S  nmi  (30  000  ft)  altitude.  The  radar  platform’s  velocity  vector  points  toward  the  top  of  the 
figure.  The  frequency  for  zero  azimuth  and  zero  dive  angle  is 
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Fig.  A2  —  Isodops  for  a  S  nmi  (30  000  ft)  altitude,  0*  dive  angle  and 
64  filters;  the  transmit  frequency  is  3  OHz;  and  the  PRF  is  SO  kHz 

2v, 

/o  —  —jp  cos  (A4g) 

For  180°  azimuth  and  zero  dive  angle  the  frequency  is 

2v 

/no  - - cos  4>.  (A4b) 

In  our  case  the  azimuth  of  the  antenna  is  90s,  v,  -  300  knots,  X  —  0.1  m,  and  4>  —  4.3s 
(corresponding  to  a  3s  grazing  angle). 

Substituting  these  values  in  Eq.  (A4a)  gives  a  frequency  at  zero  azimuth  of  3080  Hz  (-3080  at 
180s  azimuth).  The  frequency  at  90s  azimuth  and  zero  dive  angle  by  Eq.  (A3)  is  0.0  Hz. 

Also  for  our  case,  A/*  —  gives  A/*  “  781  Hz  where  N  (the  number  of  pulses  processed)  is 

.  NT 

64,  T-  2  x  nr5  s  (i.e.,  PRF  -  50  kHz). 

The  frequency  feeding  into  a  filter  is  computed  by 

fj  -  nUft  (A5) 

where 

»  -  INT(/,W) . -2,  -1,  0,  1,  2 . INT(/0), 
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and 

j  —  (n  +  1)  modulo  ( N ). 


Since  1NT  -  is  3,  there  are  seven  filters  such  that  — — —  cos  4>  <  fj  <  ~~  cos  <f>  (i.e., 

I  781  Hz  J  A  A 

there  are  seven  filters  that  contain  clutter).  In  our  case  n  —  -3,  -2.  — 1,  0,  1,  2,  3.  Filter  1  is  cen¬ 
tered  at  0.0  Hz,  filter  2  at  0  “  2,  n  -  1)  781  Hz,  filter  3  at  O'—  3,  n  —  2)  1562  Hz,  and  filter  4  at 
2343  Hz.  Filter  5  would  be  at  3124  Hz,  but  since  the  frequency  at  0°  azimuth  is  3080  Hz  (for  this 
case),  there  is  no  filter  5.  Also  by  symmetry  (in  this  case),  filters  64,  63,  and  62  are  centered  at  -781 
Hz,  -1562  Hz,  and  -2343  Hz,  respectively. 

As  shown  by  Fig.  A2  equal  steps  in  frequency  (A/d)  produce  unequal  angular  steps. 

In  the  example  presented  here,  only  seven  filters  (1,  2,  3,  4,  62,  63,  and  64)  centered  around  zero 
Doppler  contain  clutter  as  shown  by  Fig.  A2. 

Figure  A3  shows  the  two-way  circular  antenna  pattern  used.  It  has  a  1.5-m  aperture  with  uniform 
weights  on  transmit,  34-dB  truncated  Gaussian  weights  on  receive,  and  a  nominal  frequency  of  3  GHz. 


Fig.  A3  —  Two-way  circular  antenna  pattern  with  a  1.5-m  aperture,  uniform  weights  on  transmit, 
34-dB  truncated  Gaussian  weights  on  receive,  and  frequency  of  3  GHz 


To  get  the  clutter  contribution  along  an  isodop  we  must  sum  over  many  ambiguous  slant-ranges. 

To  determine  the  clutter  level  in  each  range-gate  we  sum  the  contribution  at  all  azimuth  angles 
received  along  a  constant  slant-range  cell  plus  all  ambiguous  slant-range  celts  (we  have  a  range  ambigu¬ 
ous  waveform). 

The  ambiguous  range-gates  with  equal  slant-range  increments  corresponding  to  a  50  kHz  PRF  are 
shown  in  Fig.  A4.  When  Fig.  A4  is  compared  with  Fig.  A2,  it  is  evident  that  there  are  many  different 
Doppler  frequencies  at  each  particular  slant-range. 
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Fig.  A4  —  Ambiguous  range  gales  for  a  S  nmi  (30  000  ft)  altitude 
and  PRF  of  SO  kHz 


Therefore,  to  model  the  clutter,  vertical  angles  relative  to  the  antenna  depression  angle  of  the 
beam  center  are  computed  for  each  ambiguous  slant-range  cell  starting  with  the  minimum  ambiguous 
slant-range  cell. 

Azimuth  points  relative  to  the  horizontal  pointing  angle  of  the  antenna  are  then  computed  for 
each  point  on  the  circular  antenna  pattern  (Fig.  A3).  Each  different  depression  angle  has  this  set  of 
azimuth  points. 

The  amplitude  of  each  clutter  sample  is  computed  for  each  point  on  the  antenna  pattern.  The 
amplitude  consists  of  the  gain  of  the  antenna  pattern  at  that  angle  and  the  radar  cross  section  of  the 
clutter  which  has  a  Rayleigh  distribution  in  this  example. 

Next  the  Doppler  frequency  for  each  clutter  sample  is  determined  by  Eq.  (A3).  The  total  clutter 
amplitude  at  a  given  Doppler  frequency  is  determined  by  summing  the  amplitudes  of  all  of  the  clutter 
samples  with  the  given  Doppler  frequency  in  a  selected  slant-range  cell.  This  includes  all  of  the  ambi¬ 
guous  slant-ranges  of  the  selected  slant-range  cell. 

After  correcting  for  platform  motion  by  a  TACCAR  processor  technique  [All,  the  amplitudes  for 
each  Doppler  frequency  can  be  plotted  to  construct  the  clutter  spectrum  (Fig.  2). 
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Clutter  Level 


According  to  Skolnik  (A2], 


®3dB  “ 


58.5  A 


(A6) 


for  a  circular  distribution,  where  X  is  the  transmitted  wavelength,  a  is  the  aperture  width,  and  03  dB  is 
the  hair-power  beamwidth  in  degrees. 


By  Nathanson  [A3],  the  area  of  the  clutter  is 


Ac  —  ~  sec  tan 


03  dB 
2 


(A7) 


where  R  is  the  slant-range  to  the  center  of  the  clutter  cell,  r  is  the  pulse  width,  c  is  the  propagation 
speed,  and  is  the  grazing  angle. 


The  slant-range  is  determined  by 

R  -  R,  sin  4  -  -Jr?  -  R }  cos*  <f»,  (A8) 

where 

Rp-  Re  +  h.  (A9) 

Here,  h  is  the  altitude,  Rt  is  the  earth’s  radius,  and  <f>  is  the  depression  angle  of  the  earth  point  with 
slant-range  R. 

To  compute  the  total  clutter  cross-section, 

-  Ac<t°N„  (A10) 

where  <rc  is  the  total  clutter  cross-section,  «r°  is  the  clutter  cross-section  per  unit  area,  and  Na  is  the 
number  of  slant-range  ambiguities. 


The  signal  to  clutter  ratio  is 


g| 

Ac<r°N'  ' 


(-<r°)dB  -  10  log  Ae  -  10  log  N,  +  (<r,)dB, 


where  a,  is  the  target  mean  radar  cross-section. 


(All) 


In  our  example  we  consider  an  5-band  (3  GHz)  radar  with  a  1.5-m  aperture  at  an  altitude  of 
5  nmi  (30  000  ft).  The  range  resolution  (cr/2)  is  90  m  (300  ft).  For  this  case,  W,- 102,  X- 0.10m, 
«  ■  1.5  m,  ♦  ■  3#,  and  <r,  “  13  dBsm. 


Equation  (A6)  produces  03  &  -  3.90*.  Equation  (A7)  gives 

Ac  -  0.26  nmi  (9.08  x  10s  m2). 
Equations  (A10)  and  (All)  produce 


-  67  dB. 
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Table  A1  summarizes  the  results  for  various  clutter  types  where  SCR  is  the  signal  to  clutter  ratio. 
Our  goat  in  this  research  was  a  70  dB  improvement  factor  (IF). 


Table  Al  —  Clutter  Estimation 


Grazing  Angle 

SCR  (dB) 

■taiisiMijj 

Sea 

-46 

3* 

WEM 

34 

Desert 

-40 

10° 

40 

Farmland 

-33 

10° 

-34 

47 

Woods 

-26 

10° 

-41 

54 

Cities 

-18 

10° 

-49 

62 

'Reference  A3,  pp.  23$,  263 
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